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Abstract 

A recently proposed algorithm to obtain global solutions of the double 
confluent Heun equation is applied to solve the quantum mechanical prob- 
lem of finding the energies and wave functions of a particle bound in a 
potential sum of a repulsive supersingular term, Ar~ 4 , plus an attractive 
Coulombian one, -Xr _1 . The existence of exact algebraic solutions for 
certain values of A is discussed. 

Supersingular potentials, i. e., potentials presenting at the origin a singular- 
ity of the type r~ a , with a > 2, were firstly discussed in the context of collision 
of particles [1- 6] , the main issue being the adequate definition of the S matrix. 
On the other hand, concerned with bound states, a considerable interest on re- 
pulsive singular potentials added to an attractive regular one was arisen by the 
seminal paper of Klauder [7j , where the today known as "Klauder phenomenon" 
was reported. 

As shown by Detwiler and Klauder [5], supersingular potentials cannot be 
treated by conventional WKB or perturbative methods. In view of this, a 
great diversity of approximate methods have been suggested. Extensive lists 
of articles applying those methods to supersingular potentials can be found in 
recent papers by Saad, Hall and Katatbeh [9] and by Liverts and Mandelzweig 

The validity of an approximate method can be inferred from comparison of its 
results with those obtained with an exact method. Up to now, direct numerical 



javicrammzar . es 



1 



integration of the Schrodinger equation has been the only "exact" method to 
compare with. Obviously, an algebraic exact method would be preferable. 

In order to better elucidate the effect of the supersingular potential on the 
energies of a particle bound in an attractive potential, it is much safer to choose 
for this one an exactly solvable form. For this reason, the attractive part of the 
potential has been taken, in most cases, to be a harmonic oscillator or a Coulom- 
bian well. Nevertheless, even in those simple cases, the resulting Schrodinger 
equation presents two irregular singularities, one at the origin and the other at 
infinity, and solving it requires to deal with the non-trivial issue of connection 
between their singular points. 

The problem of connecting singular points of a differential equation has been 
considered by several authors. Most of them [11-14] refer to a differential equa- 
tion for which the origin is an ordinary or a regular singular point and the 
infinity is an irregular singular one. As far as we know, only the procedure 
developed by Naundorf [15] becomes applicable also to the case of both the 
origin and the infinity being irregular singular points. In a recent paper [16j . an 
algorithm, based on a modification of the Naundorf 's procedure which improves 
notably its performance, was presented and applied to find global solutions of 
the double confluent Heun equation, a simple case of differential equation with 
two irregular singular points. The purpose of this note is to show the usefulness 
of that algorithm to obtain exact solutions of the Schrodinger equation with 
supersingular potentials and, therefore, to test the quality of the different ap- 
proximate methods proposed to deal with them. We concentrate on a specific 
potential, namely an attractive Coulomb potential to which a singular repul- 
sive term has been added, whose discrete energy levels have been studied by 
Aguilera-Navarro et al [17] by using variational methods. 

The quantum mechanical problem to be solved is that of a particle of mass 
m and angular momentum hi bound in a spherically symmetric tridimensional 
potential whose radial dependence is given by 

^>=£(4?-^f)- A ><>- z ><>- « 

Here ro represents an arbitrary length fixing the scale. The factor h 2 /2m in 
front of the parenthesis guarantees that the parameters A and Z, that represent 
the intensities of the supersingular and the Coulomb parts, are dimensionless. 
The energy of the particle will be denoted by 



2mr\ 

where E is also dimensionless. We are here interested on the bound states of 
the particle. In other words, we are looking for values of the parameter E < 
for which the Schrodinger equation admits normalizable solutions. This requires 
the wave function to be regular also at the origin and at infinity, the two singular 
points of the differential equation. 
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The changes of independent and dependent variables 

r/r Q =z and R(r) oc z~ 1/2 y(z) (2) 
turn the Schrodinger equation for the radial wave function, 



h 2 /l d ^ d 1(1 + 1)\ + y( yA R r r \ _ h 2 f p. ( -,, 
2m V r 2 dr dr r 2 ) ) 2mr? 



o 

into a particular case of the double confluent Heun equation [18) 

d 2 

D 2 y(z) + B(z)y(z) = 0, D = z—, B(z)=J2 B p z ^ ( 4 ) 

p=-2 

with parameters 

B_ 2 = -A, B_i = 0, B = - V 4 , Si = Z, B 2 = ^- (5) 

The algorithm described in Ref. [TB] deals with the function 

«;(*) = zy 2 y(z), (6) 
which turns out to be proportional to the reduced radial wave function, 

w(z) oc rR(r), (7) 
and satisfies the differential equation 

z 2 ^+(-^-l(l + l) + Zz + Ez 2 ^jw = 0. (8) 

Now it is immediate to use our method to determine the eigenvalues E and the 
(reduced) eigenfunctions w. 

The differential equation Q admits two Floquet solutions whose Laurent 
expansions are 

oo 

Wj = z v i c «>^"' 3 = ^X (9) 

n— — oo 

the coefficients c n j obeying a fourth order recurrence relation, 

—Acn +2 ,j + ((n + Vj)(n - 1 + Vj) - 1(1 + 1)) c nj - + Zcn-i,j +Ec n - 2 ,j = 0. 
The requirement of those expansions to be convergent, i. e., the condition 

oo 

|c nj | 2 <oo, (10) 

n=— oo 

fixes, for given values of the parameters A, I, Z, and E, each one of the indices 
Vj , except for addition of an arbitrary integer accompanied of a relabeling of the 
coefficients c n j . Any solution of ([8]) can then be written as a linear combination 

w — Ci wi + C2 w 2 (11) 



3 



of the two Floquet solutions. To be acceptable from the physical point of view, 
the behaviour of that combination at the singular points must be 

w(z) oc exp(— az) z M (l + 0(z -1 )) , as z — > oo, (12) 
w(z) oc exp {-0z~ x )z (l + 0(z)), as z^O, (13) 

with exponents 

a = V^E, U= — = VA. 

Although the combination pip contains two arbitrary coefficients, £i and C$,, 
the homogeneity of the differential equation reduces to only one the number of 
effective degrees of freedom of that combination. Then, for given A, I and Z, 
both conditions (fT2"]l and (fT3"]) are satisfied only for the values of E corresponding 
to bound states. 

We have applied our algorithm [16] to the determination of the energy of the 
lowest states bound in the potential (P) for the set of values of the parameter 
A considered by Aguilera-Navarro et al [17] . The Coulombian parameter Z has 
been taken, following those authors, equal to 1. The results corresponding to 
any other value of Z > 0, for instance Z, can be inferred from those shown in 
Table 1 by noting that a trivial change of scale transforms Eq. ((8]), with Z = 1, 
into 

, d 2 w f A . - ,\ 

~dz Y+ ~^ -l{l + ±) + Zz + Ez 2 \ w = 0, 

with A = Z~ 2 A and E = Z 2 E. 

Besides the values of the energy, our algorithm provides with the indices 
Vj and the coefficients c nj of the Floquet solutions and with the values of the 
coefficients £i and £2 that make the linear combination ([IT]) to be well behaved 
at the singular points, i. e., the (to be normalized) solution of the physical 
problem. In the neighbourhood of the singular points, the wave function should 
be evaluated by making use of its asymptotic expansions 

00 

w(z) r~j exp(— az) z^ a m z~ m , as z — > 00, (14) 

m=0 

00 

w(z) — exp(-0z~ r ) z^2 b m z m , as z — > 0, (15) 

m=0 

where the first coefficients ao and bo are given by our algorithm and the rest of 
them are obtained from the recurrence relations 

-2ama m = ((m - /x)(m - 1 -//)-/(/ + 1)) a ro _i - Aa m - 3 (16) 
-2 0mb m = (m(m-l)-l(l + l))b m -i + Zb m - 2 + Eb m ^ 3 , (17) 

As an example, let us consider the I — ground state in a potential of intensity 
A = 10. Our algorithm gives for the energy the value E = -0.093111277969. 
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Table 1: Energies E n i of the ground (n — 0), first excited (n = f), and second 
excited (n — 2) states of a particle of angular momentum I = 0, 1 or 2 in the 
potential given by Eq. (TTJ) with Z = 1 and different values of A. The energies 
are given in units h 2 /2mrQ. In absence (A = 0) of a supersingular term, the 
energies are the Coulombian ones, E n j = 1/4 (n + I + l) 2 . 









En t i 




Eli 


E2.1 





0.0001 


-0. 


.245429530 


-0. 


.061924800 


-0.027606954 




0.01 


-0. 


.218348228 


-0. 


.058346200 


-0.026526103 




1 


-0. 


.139037013 


-0. 


.045670342 


-0.022416709 




25 


-0. 


.077060194 


-0. 


.031979568 


-0.017358131 




100 


-0. 


.056304895 


-0. 


.026031554 


-0.014892512 


1 


0.0001 


-0. 


.062499741 


-0 


.027777693 


-0.015624935 




0.01 


-0. 


.062475231 


-0 


.027769650 


-0.015621471 




1 


-0. 


.060877715 


-0. 


.027256429 


-0.015398665 




25 


-0. 


.050470052 


-0. 


.023926758 


-0.013937057 




100 


-0. 


.042059190 


-0. 


.021077780 


-0.012640912 


2 


0.0001 


-0. 


.027777774 


-0. 


.015624998 


-0.00999998 




0.01 


-0. 


.027777435 


-0. 


.015624841 


-0.009999914 




1 


-0. 


.027743922 


-0. 


.015608960 


-0.009991373 




25 


-0. 


.027076025 


-0. 


.015299319 


-0.009825968 




100 


-0. 


.025711754 


-0. 


.014683511 


-0.009499543 
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For the indices of the Floquet solutions we obtain pure imaginary values, i>\ = 
—i>2 = 0.918988880508 i. So, besides of opposite, they turn to be complex 
conjugate to each other. In this case, since the parameters of the potential 
are real, the coefficients c nj - of the Laurent expansions © of the two Floquet 
solutions are complex conjugate, c Ui 2 = c n ,i, and Wziz) = w\{z) for real z. To 
give an idea of the magnitude of the coefficients of the Laurent expansion, we 
report in Table 2 some of those of w\ normalized in such a way that c n fl = 1. A 
physically acceptable solution results if we take a linear combination as in pip 
with 

Ci = -0.059000052486 + 0.998257979586 i, (2 = - Ci • 
The coefficients ao and bo of the asymptotic expansions (|14p and (fT5| are then 

a = -2.5451305418 i, b = -4.1825205880?. 



In the above discussed example, the indices Vj of the Floquet solutions hap- 
pen to be pure imaginary. This is not always the case. We report, in Table 
3, those indices for the most tightly bound states, with different angular mo- 
menta, in potentials corresponding to several values of A. It is interesting to 
notice that, given the indices Vj vanish for certain values of A. In these cases, 
the wave function could be written in the closed form 

w(z) — exp(— az — fiz^ 1 ) z v(z), (18) 

where v(z) represents a series of integer powers of z, let us say 

00 

The question arises if this series could in fact turn to be a polynomial for any 
values of A, that is, if w(z) is a quasi-polynomial solution, in the terminology of 
Schmidt and Wolf [18]. Comparison of (fTS)) with (fl"4j) and (fT5l) makes evident 
that /1 should be an integer, for instance p. Accordingly, the energy parameter 
would be 

Z 
V 

whereas the parameter A should be such that the coefficients £j in (fT9|) vanish 
for j > p. Obviously, the £j are related to the coefficients a m of the asymptotic 
expansion (TJJJ by 

p— 1 00 
exp(-/3z- 1 ) J2 6- j = ^ E fl ™ z ~ m ' ( 21 ) 

j—0 rn—O 

and to the b m of (fTS")) by 

»— 1 00 



E =-T~2 (2°) 



cxp(-az) £ ^ z j - 2 bm z " 1 - ( 22 ) 

j=0 m=0 
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Table 2: Real and imaginary parts of the coefficients c n ,i of the Laurent expan- 
sion, Eq. (JOj) , of one of the Floquet solutions, wi, of the differential equation {8} 
with A = 10, I = 0, Z = 1, and E = -0.093111277969. The normalization of 
wi has been arbitrarily taken such that co.i — 1. The index (see the text) turns 
to be v\ = 0.918988880508 i. The other Floquet solution, ui2, has coefficients 
Cm. 2 = c„,i and index ^2 = t^T- 



71 


3?c n) i 




3c nj i 




-20 


-0.121469106359E- 


-09 


0.131600927337E- 


-09 


-19 


-0.597147461925E- 


-09 


-0.103014061043E- 


■08 


-18 


-0.459666139538E- 


-08 


0.597145217189E- 


■08 


-17 


-0.263452992860E- 


-07 


-0.369047564493E- 


■07 


-16 


-0.136571713612E- 


-06 


0.219244950133E- 


■06 


-15 


-0.923098009638E- 


-06 


-0.104082342414E- 


■05 


-14 


-0.304091218966E- 


-05 


0.635537646613E- 


■05 


-13 


-0.250549707024E- 


-04 


-0.222398101994E- 


■04 


-12 


-0.467558678970E- 


-04 


0.140924281725E- 


03 


-11 


-0.509362942365E- 


-03 


-0.340072927974E- 


03 


-10 


-0.404150379497E- 


-03 


0.229165431769E- 


■02 


-9 


-0.740381881993E- 


-02 


-0.336931518078E- 


■02 


-8 


-0.394128245178E- 


-04 


0.257592989152E- 


■01 


-7 


-0.719278491963E- 


-01 


-0.168793007782E- 


■01 


-6 


0.392262032139E- 


-01 


0.182994779417E4 


-00 


-5 


-0.419924194846E- 


fOO 


0.866008349063E- 


■02 


-4 


0.372866041301E- 


fOO 


0.704332666002E+00 


-3 


-0.121096156077E- 


fOl 


0.468201746853E4 


-00 


-2 


0.125443087271E- 


fOl 


0.103994993070E+ 


-01 


-1 


-0.100849614575E- 


fOl 


0.137165533736E4 


-01 





0.100000000000E- 


fOl 


O.OOOOOOOOOOOOE-f 


-00 


1 


0.398351667070E- 


fOO 


0.536163740419E4 


-00 


2 


-0.196983836993E+00 


0.355835389786E- 


■01 


3 


0.264749708403E- 


-01 


-0.214447855751E- 


■01 


4 


-0.204690743398E- 


-02 


0.342012088434E- 


■02 


5 


0.953365922946E- 


-04 


-0.324605520121E- 


03 


6 


-0.191766489961E- 


-05 


0.227409233686E- 


■04 


7 


-0.102980463212E- 


-06 


-0.125761848401E- 


■05 


8 


0.131145372215E- 


-07 


0.579282461546E- 


■07 


9 


-0.819888649623E- 


-09 


-0.228010314464E- 


■08 


10 


0.382880819640E- 


-10 


0.785818186305E- 


■10 


11 


-0.147533767446E- 


-11 


-0.240418109468E- 


■11 


12 


0.491035084921E 


-13 


0.662078837429E- 


13 


13 


-0.144702107807E- 


-14 


-0.165529043300E- 


■14 


14 


0.384193021031E 


-16 


0.379059006228E- 


■16 


15 


-0.929783698355E- 


-18 


-0.799671152250E- 


■18 


16 


0.207085395184E- 


-19 


0.156351350772E- 


19 


17 


-0.427466376462E- 


-21 


-0.284451924206E- 


21 


18 


0.822932028905E- 


-23 


0.483549589235E- 


23 


19 


-0.148466844882E^ 


-m 


-0.770120556643E- 


■25 


20 


0.252142755003E- 


-26 


0.115225564234E- 


■26 



Table 3: Indices v\ of one of the Floquet solutions of Eq. ((5J) with Z = 1, the 
shown values of A and I, and the energy E of the most tightly bound state. The 
indices of the other Floquet solution are the opposite, V2 = —v\. Notice that 
V2 = v\ if v\ is pure imaginary. In the cases of complex v\, since $tv\ = 1/2, the 
ambiguity in the definition of the indices allows one to take V2 = — v\ + 1 = v\- 
This choice makes the set of coefficients of the Laurent expansion of one of the 
Floquet solutions to be the complex conjugate of the set of the other one. 



1 


A 


V 




E 





0.5 


0.487816983037 i 


-0 


.1533318066 




5 


0.825876276087 i 


-0. 


1062948831 




65 


0.5 + 0.556566003844 i 


-0 


.0622769642 


1 


0.5 


0.020895739155 


-0. 


.0615761515 




5 


0.229880095066 


-0. 


.0574288855 




65 


0.5 + 0.898656568840 i 


-0 


.0448293862 


2 


0.5 


0.000919514453 


-0. 


.0277607452 




5 


0.008595680010 


-0. 


.0276154597 




65 


0.108838373090 


-0. 


.0262524928 



The relations 

C P -i-k=a k -^2 i — ^P-i-k+s, k = 0,l,...,p-l, (23) 

s=l 

following immediately from Eq. (|21[) . allow one to write explicit expressions of 
the £j in terms of /3. Substitution of these expressions in the relation 

Q = a p -J2^rts-U (24) 

8=1 S ' 

following also from (|2ip . produces a polynomial equation for f3 whose solution 
provides the values of A = f} 2 for which the wave function is given by (fl~8]) , v(z) 
being a polynomial of degree p — 1. Table 4 shows several sets of values of A 
and I for which the energy E and the wave function w of a (ground or excited) 
state admit such particularly simple expressions. Solutions of that kind, i. e., 
quasi-polynomial solutions, are especially interesting to test the accuracy of the 
different approximate methods. The same results are obtained starting from 
Eq. (|22|) instead of (j2Tj) . The relations defining the £j are then 

& = &k-£)^-&-' fc = 0,l,...,p-l, (25) 
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Table 4: Sets of values of the parameters A ^ 0, I, and E for which Eq. {5J 
with Z = 1 admits a particularly simple solution of the form (|T8|l . u(z) being 
a polynomial. The list could be continued, the values of E being given by 
E = — l/Ap 2 , with integer p. 



E 



-1/16 
-1/36 



6 (4 
6 (A 



V7) 
V7) 



30 



-1/64 



(16/3) (lO + K{ [2(436 + i 9^4434)] 1/3 V 
(8/3) (20 - ft { (1 + i V3) [2(436 + i 9V4434)] 1/3 }) 
(8/3) ^20-K|(l-»V3) [2(436 + i9V4434)] 1/3 }) 

8(8 + 3 v / 2) 
8 (8 - 3a/2) 



8 4 



22) 



-1/100 



5 20 



310 



e 2 


+ 1080/e 


1/2 


e 2 


+ 1080/e 


1/2 


e 2 


- 1080/e 


1/2 


e 2 


- 1080/e 


1/2 



310 + 25ft I [25 (231709 + i 108V1810371)] 1/3 j) /3 



1/2 



10 (16 + V37) 
10 (16 - \/37) 
30 



(10/3) (25 + 23ft I [4 (439 + i 9V1A559)] 1/3 j 
(10/3) + [A (439 + i9\/T4559)] 1/3 | 



(10/3) (lC 1 + [3142 - 9715519] 1/3 + [3142 + 9VT551~9 



1/3 



and the equation whose fulfilment determines turns out to be 



p 




,p—s ■ 



(26) 



s=l 



Notice that this equation expresses the fact that £ p , defined as in (|2l)|) with k = p, 
vanishes. On the other hand, substitution of (|18p and (|19[) in the differential 
equation ([5]) reveals the recurrence relation 



- 2(3j^ = - 1) - 1(1 + 1) - 2a(3) + (Z - 2a(j - 1)) ^_ 2 . (27) 



Therefore, since 2a = Z/p, all £j with j > p vanish if £ p = 0. By the way, the 
last recurrence relation provides a third procedure to determine the values of 
A for which a quasi-polynomial solution exists: start with £o and £i as given 
by (|25|) . use the recurrence |27|) to write £ p in terms of /3, and then solve the 
equation £ p = 0. 

The results shown above have been obtained by using double precision FOR- 
TRAN codes. Additional digits in the results would require higher precision 
arithmetic. This is especially necessary in the case of excited states, due to the 
progressive decreasing of the energy gap between consecutive levels. 
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